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ABSTRACT 

Combining the Wide-Field Infrared Survey Explorer All- Sky Release with the 2MASS Point Source 
Catalog, we detect a nearby, moderately metal-poor stellar debris stream spanning 24° across the 
southern sky. The stream, which we designate Alpheus, is at an estimated distance of ~ 1.9 kpc. Its 
position, orientation, width, estimated metallicity, and to some extent its distance, are in approximate 
agreement with what one might expect of the leading tidal tail of the southern globular cluster NGC 
288. 

Subject headings: globular clusters: general — globular clusters: individual (NGC 288) — Galaxy: 
Structure — Galaxy: Halo 



1. INTRODUCTION 

Wide-field photometric surveys have enabled the 
identification of at le ast 15 st e llar debris streams 
in the Galactic halo (jGrillmairl I2Q1QI : iBonaca et al.l 
120 121 : iMartin et all 120131 ^ Roughly the same num- 
ber of cold s ubstructur es have been found i n ve- 
locity space (iHe lmi et al. 1999; Smit h et al.l I2QQ9I : 
iSchlaufman et an.2QQ9. : .Williams et al.. .2QTTI ). In ad- 
dition to helping us reconstruct the accretion history 
of the Galaxy, they can al so serv e as se nsitive probes 
of th e Galactic potential (iLaw. M aiewski . fc JohnstonI 
l2QQ9l : iKoposov. Rix. fc Ho gg 2010)." Globular clus- 
ter streams are particularly interesting in this regard 
as they are the coldest stel l ar substructures yet dis- 
covered (Combe s etHI 119991 : lOdenkirchen et al.l 120091 : 
IWillett et al. 20091). A large sample of such streams will 
ultimately enable us to map the distribution of Galactic 
dark matter with much greater spatial resolution than 
has been possible to date. Detailed analyses of these 
streams may also enable us to better characterize the 
dark mat ter sub-halos believed to populate the Galac- 
tic halo. ( lYoon, John ston, fc Hoggll2010l : lCarlb"er3 120091 : 
iCarlberg et al. 2012). 

In this paper we examine the 2MASS Point Source Cat- 
alog (Skrutskie et al.. 2006) and the Wide Field Infrared 
Surve y Explorer (WISE) All-Sky Release (IWright et al.l 
I2OIOI ) for evidence of nearby streams in regions of the sky 
as yet unexplored in other surveys. We briefiy describe 
our analysis in Section [21 We characterize the stream in 



Section[3]and we discuss evidence that the stream may be 
associated with the globular cluster NGC 288 in Section 
13. li We make concluding remarks Section HI 

2. DATA ANALYSIS 

Data comprising WISE bands Wl (3.4/im) and W2 
(4.6/im), along with matching 2MASS J (1.25/im) pho- 
tometry, were extracted from the WISE All- Sky Re- 
lease and anal yzed using the ma tche d filter te c hniqu e 
described by RJock osi et all (|2002D and iGrillmairl (J2009D . 
Matched-filtering has been used to detect halo streams 
to ver y low surface densities in Sloan Digital Sky Survey 
data (iGrillmair fc Johnso n 2006 ; Grillmair fc Dionatosl 
l2006i a l&-illmairl I2009L I20T1I : IBonaca et all l20T2h . 
While 2MASS data have been used to map the Sagit - 
tarius stream using M giants (jMajewski et al.l l2003f ). 
both WISE and 2MASS data are more limited in this 
application for main sequence stars due to their much 
fainter magnitudes. SDSS experience has shown that 
matched-filter techniques over large areas work well only 
if the main sequence turnoff of a population is de- 
tectable. Moreover, while WISE goes somewhat deeper 
than 2MASS J, the W1-W2 colors of stars are very nearly 
degenerate compared with the photometric uncertainties. 
Combining WISE bands Wl and W2 with 2MASS J im- 
proves the situation, providing both color discrimination 
and greater depth than is possible with 2MASS photom- 
etry alone. With a limiting magnitude of J ~ 16.5, the 
WISE and 2MASS data can probe main sequence popu- 
lations out to a distance of ^ 4 kpc. 



An additional issue with the WISE Ah-Sky data are 
the presence of significant, magnitude-dependent, pho- 
tometric biases at faint magnitudes resulting from an 
overestimate of background level^. Since these biases 
can exceed 30% at the faintest magnitudes, this has seri- 
ous consequences for filtering of sources based on color. 
Aperture magnitudes are less affected by this bias but 
they are much noisier due to aperture contamination. We 
consequently measured this bias in Wl and W2 by com- 
paring profile-fitting measurements with aperture pho- 
tometry using several million well- fit, isolated stars in the 
region 0° < R.A. < 45°, -60° < dec < -30°. We gen- 
erated a look-up table of photometric bias versus mag- 
nitude at half-magnitude intervals and interpolated this 
table to compute appropriate corrections for profile-fitted 
photometry. 

We constructed matched filters based o n the Padova 
database of theoreti c al ste llar isochrones (;Marigo et al.l 
120081 : iGirardi et all l2010f ). All stars with < 
J < 16.5 were used, and we dereddened both WISE 
and 2MASS photometry as a function of position 
on the s ky usi ng the DIR B E/IRA S dust maps of 
iSchlegeL Finkbei ner. fc DavisI ([l998h. correct ed using 
the prescription of lSchlafiv fc Finkbeinerl (J2011i) . We use 
only stars with CC flags of 0, photometric quality flags of 
"A" , "B" , or "C" , x^ < 3.0, and the ext flag = 0, to select 
on uncontaminated, unresolved objects with SNR > 5. 
For the present search, the color-magnitude distribution 
of field stars was sampled using ~ 1.2 x 10^ stars in a 
broad region extending from the Small Magellanic Cloud 
(SMC) to the south Galactic pole. The filters were simul- 
taneously applied to the entire WISE and 2MASS cat- 
alogs, and the resulting weights were combined to yield 
all-sky, filtered surface density maps. 

Figure 1 shows the filtered star count distribution in 
the region of the south Galactic pole using a filter whose 
basis is an isochrone with Z = 0.0025, an age of 12 Gyrs, 
and is shifted so as to optimally sample populations at 
a distance of 1.9 kpc. This distribution was constructed 
by combining the weights of stars in 0.2° bins based 
on their distances from the J — Wl and J — W2 main 
sequence and giant branch color-magnitude loci. 

3. DISCUSSION 

Apparent in Figure 1 is a long, narrow feature span- 
ning some 24° across the southern sky. The orientation 
of the feature does not precisely align with lines of ei- 
ther constant R.A. or dec, or constant ecliptic longitude 
or latitude, indicating that the feature is not an artifact 
of sampling or data quality discontinuities in either of 
the two surveys. The feature shows no correlation with 
reddening maps of the area (where E(B-V) is typically 
~ 0-02), nor does it correlate with the 2MASS j_msnrlO 
valuefl which estimate the depth of the 2MASS survey 
as a function of sky position. The feature is also localized 
in color-magnitude space, which would not be the case 
if the feature were due simply to variations in complete- 
ness. We conclude that the feature is not an artifact, 
but a bona fide stellar debris stream. Other curvilinear 
features are evident at various locations across the sky, 
but at much lower signal-to-noise ratios. 



^jhttp : //wise2 . ipac . caltech . edu/docs/release/allsky/expsup/se 
^[http : //www . ipac . caltech . edu/2mass/releases/allsky/doc/sec6_ 



The visible stream runs through Phoenix, Eridanus, 
Hydrus, Tucana, and becomes lost at its southern end 
near the SMC. Given the rather blue filters we are us- 
ing and the magnitude limits of the WISE and 2MASS 
surveys, we regard it as extremely unlikely that the new 
stream could be associated with the SMC. Moreover, the 
orientation of the str eam is nearly perpendicular to th e 
SMC's orbital path (jPiatek. Prvor. fc Qlszewskill2008[ ). 
We designate the new stream Alpheus, after the river in 
the Illiad. The northern end of the stream appears to end 
fairly abruptly at (RA, dec) = (27.7°, -45.0°). Over the 
interval — 69 < 5 < —45, the path of the stream is well 
fit in celestial coordinates (to within 1°) by a polynomial 
of the form: 

a = 32.116 - 0.00256 xS- 0.00225 x S^ (1) 

Color-magnitude Hess diagrams (CMDs) for Alpheus 
are shown in Figure 2. These distributions were deter- 
mined using Equation 1 to select stars within 1° of the 
centerline of the stream. Similar regions, 4° wide on ei- 
ther side of the stream and laterally offset by 5°, were 
used to sample the field star population. Scaling the 
latter to the former by area, the CMDs in Figure 2 are 
subtractions of the two. While not strong, a turn-off and 
upper main sequence are clearly visible. For J - Wl, 
the sequence is matched reasonably well using a Padova 
isochrone having [Fe/H] = -1.0, though the uncertainty is 
quite large because the isochrones for different metallic- 
it ies are rather closely spaced at these wavelengths. The 
J - W2 CMD is less convincing, suggesting a somewhat 
bluer sequence. Based on the stronger J - Wl sequence. 
Figure 2 suggests an uncertainty on the order of at least 
0.5 dex. Deeper follow-up surveys in visual bands and 
spectroscopy of individual stars will be required to bet- 
ter constrain the metallicity of the stream. 

A number of curvilinear features are evident in Figure 
1, and one might ask how Alpheus compares in terms 
of significance. We can make use of the "T-statistic" of 
Grillmair (2009), which measures the median contrast 
along its length between a putative stream and the sur- 
rounding field. The T-statistic for Alpheus, comparing 
with the field extending the length of the stream and 
30° to the north and 13° to the south in ecliptic coordi- 
nates, is shown in Figure 3. The stream is detected at 
the 19(7 level, and the lateral full- width- at-half-maximum 
of the stream is 3.2°. At a distance of 1.9 kpc (see be- 
low) this corresponds to a spatial extent of ^ 100 pc. 
This is similar to the widths of other streams that are 
known o r believed to have been g enerated by globular 
clusters ("Odenkirchen et al.' 2003; 'Grillmair fc John sonl 
2006; Grillmair fc Dionatos 2006b; Grillmair 2009, 2011^ 
suggesting that the progenitor of Alpheus was itself a 
globular cluster. The next-most convincing overdensity 
is detected at < lOa, and we will consider this and other 
features in a forthcoming paper. 

Following iGrillmair fc DionatosI ()2006b[ ) we shift our 
filter brightward and faintward to estimate the stream's 
distance. Combining results for J - Wl and J - W2, we 
find that the signal-to-noise ratio of the northern half of 
the stream peaks at a distance modulus of 11. 5 ±0.4 mag, 
while the southern half of the stream peaks at 11.0 ± 0.5 
cfiiggg..h^h(3 uncertainties are quite large due to a combi- 
^n^;"?;^]! of the small number of stars, large photometric 



errors at the faint limit, and the relatively narrow range 
of stellar colors at these wavelengths. This would put 
the northern half of the stream at a mean sun-centric 
distance of 2.0 ± 0.7 kpc, while the southern half is at 
1.6 ± 0.8 kpc. We note that a simple ratio of filtered 
stream to background signal peaks at 2.3 kpc in the 
southern half and 3.9 kpc in the northern half. How- 
ever, the stream is not nearly as apparent or continuous 
using such a filter, nor does such a filter match the CMD 
in Figure 2. 

Summing the color-selected, background-subtracted, 
but otherwise unweighted star counts over a width of 
3.2° we find the number of stars in the stream to Wl 
~ 16.5 to be 310 ± 90. The quoted uncertainty is largely 
due to the shot noise contributed by the several thousand 
field stars with similar colors lying along our line of sight 
to the stream. The average surface density is 4 ± 1 stars 
deg~^, with peaks of ~ 10 stars deg~^. 

3.1. On the Possible Association of Alpheus with the 
Globular Cluster NGC , 



Based on its relatively low concentration and fairly 
eccentric, retrograde orbit, NGC 2 88 has a relatively 
hig:h predicted destru ction rate (iGnedin fc Ostrikerl 
19971). and both iGrillmair et all (jl995f ) and 
Leon. Mevlan. fc CombesI ()2000f ) found a significant 



population of extratidal stars in the vicinity of the 
cluster. Figure 4 shows projections of three possible 
orbits of the cluster, in tegrated using the Galactic 
model of lAllen fc Santillan (J991). All integrations use 
a cluster di stance of 8.3 kpc and a radial velocity of 
-47 kms~^ (jHarrisI Il996f ). The three different orbits 
are based on proper motions of jUa c osS = 4.68 ±0.2 
mas/yr, fis = -5.25 ± 0.2 mas/yr (Guo I1995S) ('a'), 
/ JQ, COS 6 = 4.64 ± 0.4 mas/yr, /j^s = —6.00 ± 0.4 mas/yr 
(|Dinescu et al.l I1997D ('b'), and Ur v cos 6 ^ = 3.8 =b 0.3 
mas/yr, jj^s = —8.1 ± 0.3 mas/yr (jPlatais et al.l [l998f ) 
(V). 

Two of the orbit integrations ('a' and 'b') lie nearly 
parallel to the new stream, offset to the east by ~ 2 — 3°. 
Orbit 'c' is also nearly parallel to the stream, but is offset 
~ 7° to the west. Tracing around the integrated orbits, 
the visible portion of the stream begins at ~ 4 kpc from 
the cluster and ends at ~ 6.3 kpc. Over the 24° extent 
of the stream, the leading arm of orbit 'a' predicts a he- 
liocentric distance ranging from 3.8 kpc at the southern 
end to 5.1 kpc at the northern end. Similarly, orbit 'b' 
ranges from 5 kpc to 6.2 kpc, and orbit 'c' ranges from 
7.2 kpc to 7.5 kpc. The near coincidences in the position, 
orientation, and distance gradient of the stream, partic- 
ularly with orbit 'a', suggest a possible association with 
NGC 288. 

Given the retrograde motion of NGC 288, the stream 
would constitute the leading arm of the NGC 288's tidal 
tails. We would expect these stars to have dropped in- 
wards towards the Galactic center before being carried 
forward along their new orbits. That the stream lies to 
the inside of putative orbits 'a' and 'b' is consistent with 
this picture. An offset of 2.5° at a distance of 1.9 kpc 
corresponds to ^ 80 pc in projection. Though this offset 
is somewhat la rger than the current tidal radius of NGC 
288 (~ 31 pc, iHarrisI ()1996f )). we expect escaped stars 
to have residual velocities that will ca rry them some- 
what farther from the cluster. Indeed, IGrillmair et al.l 



(|1995f ) and iLeon. Mevlan. fc CombesI (j2000f ) found ex- 
tratidal stars extending to at least 150 pc in the imme- 
diate vicinity of NGC 288. 

Our metallicity estimate for stars in the stream is also 
consistent with a measu red value of [Fe/H] = -1.24 (Z 
^ 0.0015) for NGC 288 (|Harrisl| 19961 ). As noted earlier, 
while the stream appears strongest using isochrones with 
Z = 0.0025 (Fe/H ^ -1.0), due to the smah number of 
stars and the narrow range of colors at these wavelengths, 
the stream is only marginally less apparent (S/N ^ 17.5) 
using an isochrone with Z = 0.0015. 2MASS and WISE 
colors are not very sensitive to metallicity, and given 
the additional uncertainty contributed by our substan- 
tial corrections to the WISE profile-fit photometry, we 
consider the agreement to be quite good. 

If we adopt a proper motion for NGC 288 of /i^ cos 6 = 
4.75 mas/yr, fis = —4.60 mas/yr, we arrive at orbit 'd' 
in Figure 4. This value of the proper motion was ar- 
rived at by both forcing the orbit to pass through NGC 
288, and to pass 2.5° to the east along the length of 
Alpheus. The heliocentric distance of this orbit in the 
range —65° < S < —45° is reduced to between 2.3 and 
3.1 kpc which, to within the uncertainties, is in agree- 
ment with our distance estimate for the stream. The 
adopted value of jUa cos S is in good agreement with the 
measured values, while the value of /i^ is 1.2cr from the 
mean of the three measurements. Combining the mea- 
surement uncertainties with inaccuracies in our adopted 
model for the Galactic potential, we consider the associa- 
tion of the stream with NGC 288 plausible, if not proved. 
Validation of our hypothesis will require radial velocity 
measurements for a significant number of stream stars. 

That the stream does not visibly connect with NGC 
288 in Figure 4 may simply be due to the limiting mag- 
nitude of the WISE and 2MASS surveys. At the north- 
ernmost point of the visible stream, the NGC 288 orbit 
model predicts that the heliocentric distances rapidly in- 
crease to 8.3 kpc at the position of NGC 288. As we 
use our filters to probe beyond 3-5 kpc, the turn-off and 
subgiant branch fall below our limiting magnitude and 
we are forced to rely wholly on the RGB for our signal. 
However, the RGB is much more sparsely populated and 
more highly contaminated by foreground stars than is 
the turn-off region. We detect NGC 288 itself only by 
virtue of the relatively large number of RGB stars in the 
cluster itself. 

It is also possible that the stream has been par- 
tially disrupted by a dark matter subhalo or other per- 
turber of significant mass and/or favorable trajectory. 
However, examinations of the Pal 5 and GD-1 st reams 
CCarlberg et al. | [20T1 ICarlberg fc Grillmaid[2QT3) have 
not identified gaps anywhere near 4 kpc in length. We 
consequently consider this an unlikely possibility, but one 
which will need to be tested with deeper observations. 

4. CONCLUSION 

Using the 2MASS and WISE All-Sky Catalogs, we 
identify a new, 24°-long, moderately metal-poor stellar 
debris stream in the southern hemisphere which we des- 
ignate Alpheus. At a distance of ^ 1.9 kpc, the stream 
is much nearer than other spatially mapped streams that 
have been found to date. 

Based on near coincidences in position, orientation, 
distance, distance gradient, and metallicity, we suggest 



that the stream may be the leading tidal arm of the 
globular cluster NGC 288. Testing this hypothesis will 
require both deeper imaging in the visible bands and 
spectroscopy. Deep er imaging by the SkyMapper survey 
(jKeller et al.l 120071 ) should enable us to reach well down 
the main sequence to better characterize the stream and 
extend its path across the sky. Based on existing proper 
motion measurements for NGC 288 as well as our at- 
tempt to estimate an orbit that would connect Alpheus 
with the cluster, we would predict radial velocities in the 
range —66 < Vr < —36 km s~^ at the northern end of the 
stream, and 23 < Vr < 36 km s~^ at the southern end. If 
indeed Alpheus is associated with NGC 288, these mea- 
surements will also enable more comprehensive modeling 
of the orbit of NGC 288 and the shape of the Galactic 
potential towards the south Galactic pole. 
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Fig. 1. — High contrast, filtered surface density map of WISE/2MASS stars in the south Galactic cap. The upper panel shows the map 
with a log stretch, smoothed with a Gaussian kernel of width 0.4°. Lighter areas indicating higher surface densities. The map is the result 
of a filter based on a Padova isochrone with [Fe/H] = -1.0, an age of 12 Gyr, and shifted to a distance of 1.9 kpc. I n the lower panel, 
the map has been background-subtracted using a window-smoothed version of itself (unsharp-masking, e.g. ILee et aL (2009)) to remove 
power at low frequencies, reduce the dynamic range, and make various features and landmarks more apparent (e.g. the LMC and SMC). 
The map is smoothed as above and is shown with a linear stretch. WISE-specific features resulting from moon-avoidance maneuvers and 
momentum dumps (repeated firing of the spacecraft torque rods at nearly the same ecliptic latitude) are indicated. The curved structure 
between the LMC and SMC is part of a ring-like feature surrounded the south ecliptic pole and is an artifact of image "decimation" in the 
WISE photometry pipeline, required in high-coverage areas due to computational limitations. 
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Fig. 2. — Hess diagrams of stars lying within 1° of the centerhne of Alpheus. The sohd hne shows a Padova isochrone with [Fe/H] =-1.0, 
age 12 Gyrs, and shifted to a distance of 1.9 kpc. The dashed lines show similar isochrones for [Fe/H] = - 2.4 and solar metallicity. Lighter 
areas indicate higher surface densities. 
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Fig. 3. — The lateral profile of the stream, as measured using the T statistic (arbitrary units). The profile measures the median signal- 
to-noise ratio in five, 5° segments of Alpheus, where the noise is measured in fields extending 13° and 30° to the south and north, 
respectively. 




Fig. 4. — Projections of the orbit of NGC 288 using three different proper motions measurements (see text), 
both forcing the orbit to pass through NGC 288, and to pass 2.5° to the east along the length of Alpheus. 



Orbit 'd' was generated by 



